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Background: Deuteron induced reactions serve as surrogates for neutron capture into compound states. Al¬ 
though these reactions are of great applicability, no theoretical efforts have been invested in this direction over 
the last decade. Purpose: The goal of this work is to establish on firm grounds a theory for deuteron in¬ 
duced neutron capture reactions. This includes formulating elastic and inelastic breakup in a consistent manner. 
Method: We describe this process both in post and prior form distorted wave Born approximation following 
previous works and discuss the differences in the formulation. While the convergence issues arising in the post 
formulation can be overcome in the prior formulation, in this case one still needs to take into account additional 
terms due to non-orthogonality. Results: We apply our method to the 53 Nb(d.p)X at E d =15 MeV and 25 MeV, 
and are able to obtain a good description of the data. We look at the various partial wave contributions, as well 
as elastic versus inelastic contributions. We also connect our formulation with transfer to neutron bound states. 
Conclusions: Our calculations demonstrate that the non-orthogonality term arising in the prior formulation is 
significant and is at the heart of the long-standing controversy between the post and the prior formulations of 
the theory. We also show that the cross sections for these reactions are angular momentum dependent and there¬ 
fore the commonly used Weisskopf limit is inadequate. Finally we make important predictions for the relative 
contributions of elastic breakup and non-elastic breakup, and call for elastic breakup measurements to further 
constrain our model. 
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I. INTRODUCTION 

Neutron capture reactions A(n,y)B are very important in 
astrophysics, for the production of heavy elements, but are 
equally relevant in stewardship science, since it is through 
neutron capture that fission is induced and chain reactions be¬ 
gin. Most often, for low energy neutrons, the capture pro¬ 
ceeds through continuum states, forming a compound nucleus 
(A+n—> B*) that then decays, either through gamma emission 
to the ground state, or through other particle channels. The 
direct measurement of neutron capture is challenging, partic¬ 
ularly because most of the targets of interest have short half- 
lives and neutrons cannot be made into targets. A proposed 
alternative it to use deuterons as surrogates [fl. In this indi¬ 
rect method, the proton inside the deuteron behaves mostly as 
a spectator, the neutron inside the deuteron is delivered to the 
target surface and gets absorbed by the target A(d,p)B*. Since 
this is a compound nucleus reaction, the final compound nu¬ 
cleus decays in the same way as after A+n—>B*. An example 
of a recent measurement applying the surrogate method is the 
171.173 Yb(d.py) experiment performed at LBNL 0. The mea¬ 
surement was performed on a nucleus for which neutron cap¬ 
ture cross sections were already available. The neutron cap¬ 
ture cross-sections extracted in the deuteron induced reactions 
were in fair agreement with those measured directly. 

While the exclusive process A(d,pn)A (with A left in its 
ground state) is clearly identified as elastic breakup, the rest of 
the cross section arising from the inclusive process A(d,p)X, 
where only the proton is measured in the final state, is harder 
to name because it encompasses many different processes. 
Some refer to this component of the cross section as inelas¬ 
tic breakup, breakup-fusion or partial fusion. In this work, we 


will always use the term non-elastic breakup. 

Although there are important applications of the surrogate 
method for neutron capture, no theoretical development has 
taken place in the last two decades to establish the method 
on firm grounds. In terms of direct reaction mechanisms, this 
process can be seen as inelastic breakup followed by fusion, 
or neutron transfer to the continuum. Significant theoretical 
effort took place in the eighties with the main idea being first 
introduced by Kerman and Me Voy 0 , with the works of Uda- 
gawa and Tamura (4) and Austern and Vincent 0 appearing 
shortly after. In |4[ the authors assume the process A(d,p)B* 
is a two step process, first breakup of the deuteron followed 
by fusion of the neutron. They describe this in distorted wave 
Born approximation in prior form and make a number of ad¬ 
ditional approximations (we will denote this theory as UT). 
In 0, the starting point is the post-form distorted wave Born 
approximation, and the authors assume that the target gets ex¬ 
cited only by the neutron-target interaction (we will denote 
this theory as AV). One difficulty in this method is the con¬ 
vergence of the matrix element. In both 0 and 0, all (d,p) 
transfer cross sections to the excited states are summed with¬ 
out explicitly introducing the properties of these states. This is 
of course a key aspect in describing inclusive processes. As it 
turned out, cross sections obtained with the post (AV) theory 
did not agree with those obtained with the prior (UT) theory. 
This generated a heated controversy that lasted a decade and 
was never fully resolved. 

A detailed analysis of the UT theory is presented in 0, 
where several approximations are considered including the 
zero range approximation and the surface approximation. The 
authors of 00 argue that the theory of AV includes unphys¬ 
ical components which should be corrected by inclusion of a 


non-orthogonality term. On the other hand, Ichimura et al. 
sa in their detailed examination, conclude that certain im¬ 
plicit approximations made in the optical reductions by UT 
are at the heart of the disagreement. Although several groups 
revisited the matter later EE1, establishing a relationship 
between the various theories, the controversy on the relevance 
of the non-orthogonality term was never resolved. Nowadays, 
many of the approximations made in the early eighties have 
become unnecessary, and thus it makes sense to revisit the is¬ 
sue. 


Two recent works have applied the AV theory to study a 
variety of reactions caul. In Lei and Moro EH these re¬ 
actions include deuteron breakup on >s Ni at intermediate en¬ 
ergies, deuteron breakup on 93 Nb at lower energies and 6 Li 
elastic scattering on 209 Bi around the Coulomb barrier. In or¬ 
der to deal with the convergence of the amplitude, the authors 
in El construct continuum bins corresponding to a square in- 
tegrable wave packet obtained averaging the scattering states 
over energy. These studies E3tH show that, overall, the AV 
theory provides a good description of the processes consid¬ 
ered. 


Given the recent experimental interest in the surrogate 
method, it is critical to develop a theory that is practical and 
reliable. Our overall goal with this work is exactly to estab¬ 
lish such a theory, which can then serve as a starting point 
to make further improvements in the future. In this work we 
will present our own derivations of the elastic and inelastic 
breakup amplitudes in both the post (AV) and the prior (UT) 
forms, within the distorted wave Born approximation (section 
II). We will show that these two theories are indeed equiva¬ 
lent if no further approximations are made, although the post 
formalism introduces numerical difficulties which are avoid¬ 
able by the prior formalism. We will also show that a non¬ 
orthogonality term and additional cross-terms naturally arise 
in the UT theory. This is a consequence of the fact that when 
writing the amplitude in the prior form, no easy separation be¬ 
tween the breakup process and the excitation process is pos¬ 
sible. In other words, in the prior formalism, breakup and 
excitation are entangled. These days, however, their computa¬ 
tion poses no difficulty. Following the work of El . in section 
III, we will use the same framework to study neutron trans¬ 
fer to bound states. This establishes an important connection 
between scattering and bound states, and provides a stringent 
test for the theory. We apply the method to surrogate reac¬ 
tions of deuterons on 93 Nb for beam energies in the range of 
10-20 MeV to compare with data, as done by El (section V). 
We analyze our results and dissect the various contributions, 
in terms of angular momentum, as well as the relative mag¬ 
nitude of elastic versus non-elastic cross sections. Finally, in 
Section VI, we draw our conclusions and provide an outlook 
into further possible developments. 



FIG. 1: Definition of the coordinates used in our formulations. 

II. THEORETICAL FORMULATION 
A. General formalism in the post representation (AV) 

Let us consider the reaction A(d,p)B* which includes elas¬ 
tic breakup and any other inelastic processes. In this section, 
we follow closely the work by Austern et al. {5} and derive the 
post-form amplitude for the process. We will adopt a spectator 
approximation for the proton, which means the proton-target 
interaction will not explicitly excite the target A. We will thus 
start by substituting the proton-target interaction V Ap (r Ap ,^ A ) 
by an optical potential U Ap (r Ap ). In addition, for the purpose 
of our derivation, we have considered A to be infinitely mas¬ 
sive. However, in the actual numerical applications, the recoil 
of the nucleus A is fully taken into account. 

The three-body Hamiltonian for the problem is 

ll — K n + K p + h A (^ A ) + V pn (r pn ) + V/\ n (T/\n , 4'a ) + U Ap (r Ap ), 

(1) 

where K n and K p are the kinetic energy operators acting on the 
neutron and proton coordinates respectively. We now consider 
the model wavefunction 

y =Xi<t>d ( l>A+Xf ( l >C B’ ( 2 ) 

where fid is the deuteron eigenfunction, 4>a is the ground state 
of the target nucleus, and (f> c B represents the c’th eigenstate of 
the final compound nucleus B. These wave functions satisfy 

hd4*d — (Knp Vnpj fid — Edfidi 

h A <p A = S A fi A . (3) 

httpB ~ (K-n + ll A + VAn) fip ~ B' 

Here, K np is the kinetic energy of the neutron-proton motion, 
and h A is the internal Hamiltonian of the target nucleus A. In 
Fig. B we define the coordinates that we use in our formula¬ 
tion. These will become useful throughout this section. 

Considering only first order in the couplings between the 
incoming deuteron channel i and the final proton channel /, 
the coupled equations for the unknowns^, and xf simplify to 


2 


the distorted wave Bom approximation (DWBA) differential 
equations, and can be written in either prior and post form 

03, 


(E t - Hi) Xi = 0 


(4) 


(g — \ — / prioA 1 !) rift A j Xl 

\ f- f)Xf- | (^VpoJfa^Xi + (Hf - Ef ) (f B \<p d <p A )xi, 

where the initial and final Hamiltonians are 


H, = K d + U Ad 
H f = K p + U Ap , 

with the prior and post operators 

Vprior — VAn U A p ~~ U A d, 

Vpost — Vup. 


(5) 


( 6 ) 


and the energies E, = E - s d - s A and Ef = E - s c B . For a 
target A with a big but finite mass, the post form will include 
a negligible remnant term U Ap - Ub p , Ub p being the optical 
potential between the proton and the nucleus B. In the ap¬ 
proximation in which the proton is treated as a spectator, it is 
important to note that, in any case, this term does not depend 
upon the intrinsic coordinates of the target A. 

According to eq. (J4|, an auxiliary wavefunction Xf defined 
by 


X f =Xf + (4> C B\<f>d<f>A)Xi, 

is such that, in the post representation, 

(Hf - Efj X f = {(p c B \V post \(/) d (f> A )xi- 


(7) 


( 8 ) 


The non-orthogonality term yp c B \4> d 4>A) vanishes for large val¬ 
ues of the proton coordinate, so Xf and Xf are asymptotically 
identical. Therefore, the T -matrix Tdc for the detection of the 
proton after the population of the state </> c B of the residual nu¬ 
cleus has two equivalent expressions: 

Tdc = {x ( f ) <f> C B\Vprior\<Pcl<f>AXi) = (T/ _) <Pb\ V postlMAX•) > (9) 

and the celebrated post-prior equivalence in DWBA is veri¬ 
fied. The proton distorted wave^^T 1 satisfies the equation: 


(E, - = 0. 


( 10 ) 


The inclusive cross section is summed over all B channels 
of a given energy. The inclusive breakup cross section for the 
process d + A —» p + B is thus 


d 2 cr 

dCl p dE p 


2n 

Tw d 


Pp(Ep) £ |(*/VbI V I Max )f S(E-E p -s c b ) 


—p p (E p ) 3 ( <t> d <pAXi\V^ \x\ ') G b (x ( f VI ^ d<pAXi) , 


( 11 ) 


where the round brackets indicate that we only integrate over 
the proton coordinate. Here, V stands for either V pr i or or V l 


Pp(E p ) 


nipkp 


post 9 


( 12 ) 


is the proton level density, and E p is the kinetic energy of the 
detected proton. 

The energy-conserving delta function 6{E — E p — s B ) can 
be written as the imaginary part of the spectral expansion 


G b = lim ) 


\<P C b) {<P C b 


e—>6 E — E n 


■ s B +16 


(13) 


This Green’s function, in operator form, is simply 
G b = 1 


E - E p - h B + ie 


(14) 


Since V np is not contained in h B , Eq.( 14 1 can be interpreted 


as a propagator in the breakup channel, i.e., when the deuteron 
bound state is absent. Furthermore, if we choose the post rep¬ 
resentation of Eq. GD- we have the advantage that the poten¬ 
tial V = V post = V pn (r p „) does not depend on £ 4 , so we can 
directly project the Green’s function onto the target ground 
state, the so-called optical reduction. 


G° b p ‘ - (<Pa\G b \<Pa) - 


1 


E E p G a K n U An {r An) "t l€ 


(15) 


The exact form of U An obtained from the above projection op¬ 
eration is the Feshbach optical potential for the n -A system in 
the ground state of A, and includes the effects of all possible 
excitations of A due to the interaction with the neutron. In the 
present paper we use, instead, a global optical potential ob¬ 
tained from a systematic fitting to elastic scattering data ED. 
For both positive and negative energies, this optical potential 
U A „ = V An + iW An (with V A „, W An real) should be interpreted 
as an energy-averaged approximation to the exact Feshbach 
potential. 

The inclusive cross section has thus the post-form expres¬ 
sion 


d 2 cr 

dElpdEp TiVd 


= -j^Pp(Ep)3 ( <PdXi I Vpn | X ( f } ) G° b p ' (xf ’| v pn | (pdXi 

(16) 


Because the real potential V np (r np ) cannot excite the nucleus 
A, this formalism suggests a two-step mechanism: V np (r np ) 
first breaks up the deuteron and then G° B propagates the sys¬ 
tem in the breakup channel, eventually leading to the absorp¬ 
tion of the neutron in the complex field U An . 

It is useful to extract from Eq. © the contributions of 
breakup without the excitation of the target A (elastic breakup, 
EB), and non-elastic breakup (NEB) where the target no 
longer remains in its ground state. To this purpose, we trans- 
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form 


B. General formalism in the prior representation (UT) 


G'f =G 0 (1 +G° H pt U An ) 

= (1 + G° B p,f Ul) Go (1 + Gf U An ) - Gffl^T. 

G°f =(1 + G^Ul)Gl(l + G° B pt U An )-Gfu An G° B pt , 

(17) 


where Go is the propagator in free space. We now define 
the scattering waves bf„(r„;£„)) subject to the optical poten¬ 
tial U An \ 

\Xn(r A n, k „)) = (1 + GfU An ) Uro (r A „-, k n )), (18) 

where \xo(r An ',k n )) is a free plane wave with momentum k n . 
From these expressions, we obtain 

3 Gf =(1 + Gful) 3G 0 (1 + G° B p, U An ) + cfw An G\f 

k 2 \ 

— ^ Am kn )I 

2m n ) 

+ G°fw An G° pt . (19) 

The first term of Eq. 0 - when used in Eq.(fT 6 |), gives rise 
to DWBA elastic breakup. For this we have transition matrix 
elements 

T E b( k,„ k p ) = V pn | (f>dXi) ■ (20) 

The elastic-breakup differential cross section for the protons 
is given by integrating over the neutron angles: 



Computing the neutron wavefunction Eq.(|24|> presents well 
known numerical difficulties. As the deuteron wavefunction 
(j>d and the post potential V pn have the same argument, the 
source term ( |23) i oscillates indefinitely as a function of the 
neutron coordinate, and the integral in Eq.(p5| converges ex¬ 
ceedingly slowly. There are methods to deal with this numer¬ 
ical problem. Alternatively we follow Udagawa and Tamura 
(UT) and revert to the prior representation. Let us first write 
the inclusive cross section in the prior representation, 

d 2 cr 2 

dSlpdEp tiVd^ P ^ P ^ 

X 3 ( MaXi\ Verier I*/'’) G B (^'| V prior \ <pd<f>AXi) > 

( 26 ) 

where V pr i or = Va„ + Ua p - Ud is the prior interaction potential 
(see Eq.(|6])). Now, however, the optical reduction defined by 
Eq. ( p3] > cannot be readily applied, because V pr j or acts on the 
intrinsic coordinates g A of the nucleus A and does not com¬ 
mute with 4>a- 

We will thus split V prior in a term that can induce breakup 
(and commutes with <f> A ) and a term that can excite 4>a (and 
commutes with X^X 

Vprior = (G b 1 +V prior )-G- 1 

= ( U A p - U Ad - h A - K n - E p + E)~ G b x . (27) 


d 2 cr 


dEpdElp 


EB 


2n 

Twd 


Pp(E p )p n (E„) J' \T EB (kn, kp)| 2 dQ. Bn , 


( 21 ) 


where, in conjunction with the proton density introduced in 
Eq.{T 2 |. we use the neutron density of states 

/z? \ m nkn 

Pn( n) ~ (22) 

The second term of Eq.(p~9} gives rise to the non-elastic 
breakup cross section. Defining the post form of the source 
term, 

S post = (*/’| V pn | <f>dXi) , ( 23 ) 

and the neutron wavefunction 

I Vn OSt )=G° pt S P osu ( 24 ) 

the non-elastic breakup differential cross section can then be 
written as 

NEB 9 

= - J-Pp( e p) (<Ar'| W An \ilr p r'} ■ ( 25 ) 

Eq. ( fl9| ) suggests that the elastic and non-elastic breakup 
have a common origin. In Appendix A, we show that flux 
conservation laws give rise to separate elastic and non-elastic 
terms in a consistent manner. 


d z cr 

dClpdEp 


Note that, when doing the averaging over the states of the tar¬ 
get, some terms drop out. Introducing A U - U Ap - U A d , we 
can write: 

{H\Vlri„ r \Xf)G B {x^\V prio rW A ) 

= (A rf - K n + E„) ffi) G opt (xfl (A u - Kn + E„) 

+ 1 X?)(Kn ~ Va„ - E n )(x^\ + (A f/ f + AU)\x^Wf ] \ 

(28) 


where E„ = E - E p - s A and V An = (<p A \V An \(l> A ). When taking 
the imaginary part in Eq.( 26 1 , the last term in Eq.( 28 i gives no 
contribution since it is a real operator. If we substitute E„ - K„ 
with UAn- we have 


3 ( M V lr ior 1 X ( f ’ )Gb (x ( f \Vprior I 4> A ) = 

= (Ac/ + + U An ) kf^GT^l (A U + U An ) 


+ \x ( f\U A n - V An )(x ( f\ (29) 


It is convenient to rewrite Eq|29]in the following form: 

3 ( MVU [ X?)G B [X^Wprior I <P A ) 

= (A U + U A n ) f Wf^GTtXf) (A U + U An ) 
+2iW An \x ( f ) )G° B p, (Xf ) \ (A U + UAn) 
+\X ( f ) )(U A „ ~ V A n)ix ( j\ (30) 
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The substitution of E n - K n with U An might seem dubious, 
because, even if it is clear from the first line of Eq. ( |TT) that we 
are dealing with on-shell quantities, the second line seems to 
formally include off-shell terms. But, if we consider E n -K n = 


G 0 !' + U An , we can see that the off-shell term G°f' 


gives 


a real contribution to Eq.(|28), and vanishes when taking the 
imaginary part in Eq.d 


It is now convenient to define the prior-form source term 


S prior = (XfYU + U An \(/> d Xi) • ( 31 ) 


Note that the operator AU + U An differs from the prior inter¬ 
action defined in Eq.(|6), namely the complex optical potential 
U An is used instead of the real interaction V An that causes core 
excitation. The prior-form neutron wave function is 


,i,P rior - n opt K 

Yn ~ ^B ^ pri° r ' 


(32) 


We also need to introduce the non-orthogonality function: 


K 


= W1 


<f>dXi) ■ 


(33) 


This last expression defines the neutron “source” function ob¬ 
tained by Hussein and McVoy (HM) using somewhat different 
approximations QoJ. Now, inserting ( [30] ) in ( |26) , we can ar¬ 
rive at: 


d~cr 

dQ.pdE p 


^(E p )[3(s pn o I \G° B pt \S pnor ) 


+ 2% (<C*| W All G° B pr |S pr i or ) + ( irf M \ W An \itf M )\, ( 34 ) 

We can apply the identity ( p~9] > to remove from the first term 
of (34 1 , the elastic breakup contribution. The remaining terms 


represent the total non-elastic breakup cross section: 


_Sa_ 

dElpdEp 

+ 2%(^ 


NEB 


= -^Pp(E p )[d (i/, prior \W An \M 


prior\ 


) 


HM \ W An \K ri ° r ) + (<A n M \W A „ \tf M )\ • 


(35) 


The first term corresponds to elastic breakup followed by cap¬ 
ture, while the third one contains all other processes involving 
the n + A system. As we will show all three terms are, in 
general, important, and have to be simultaneously taken into 
account. For completeness, we present in Appendix B the par¬ 
tial wave decomposition of these results. 

We have just shown that the non-orthogonality term in the 
above expression arises when disentangling the elastic and 
non-elastic breakup contributions from V pr j or . On the other 
hand, an identical non-orthogonality function appears for dif¬ 
ferent reasons in the derivation of the standard DWBA equa¬ 
tions in the post form, but has been dropped because of the 
equivalence of the final (proton channel) distorted waves Xf 
and Xf in the asymptotic region (see eq. 0). Nonetheless, 
in Eq.( 25) the proton distorted wave is certainly needed for 
small values of the proton coordinate, so the question may 
arise whether the latter non-orthogonality function should be 
kept after all for its derivation. 


To see that such a term is actually not needed, let us check 
the consequences of including it. By examining eq. 0, it can 
be seen that this is equivalent to making the replacement: 

V -> V p „ + (. K p + U Ap + h B - E) (36) 

in eq. ( fTT) . We obtain from the post-form matrix element in 
Eq.lfTT}: 

(<Pd<p A Xi\Vpn + ( K p + u\ p + h B - e) \x ( f ^)g b 
x (Xf J l (Vpn + (Kp + U Ap + h B - £)) 1 4>d(f> A Xi) = 

(M [S* p0st + ^ M * (h B + E p -E)\G B 
x [spost + (hB + Ep - e)] | (f ) A ), 

where we have used eqs. ( [TO) and ( |24[ >. Taking into account 
that (h B + E p - E^G b - 1, we get 

S; ost G opt S post +2%(s post ^ M ) 

+ S* post (li B + E p - e) 5 p OS t | (f> A ) ■ (37) 


When taking the imaginary part, only the first term of the 
above expression survives, and we obtain again Eq.( 16 1 . We 
conclude that the post source term defined in eq. (23 1 should 
not include a non-orthogonality contribution. 


C. Transfer to bound states 


The sort of breakup processes we are considering can be 
thought of as transfer to the continuum. One would like to 
have one framework to describe both transfer to bound states 
and continuum states. The Green’s function formalism al¬ 
lows for this connection. This idea was first proposed in 03) . 
We write the partial wave coefficients of the Green’s function 
(r An , r ' An ) defined in Eq.( 151, as 


r °pt 

G B 


Gi(r An , r’) = 


fl(k n , r An< )gi{k „, r An> ) 


k„r An r’ 


(38) 


An 


where k„ = y/2m„e/ti, and f t (k„, r An ) ( gi(k n , r An )) is the regular 
(irregular) solution of the homogeneous equation 


n 2 


2m » dr ln 


U A n(r A n) 


n 2 i(i+ 1) 
2m An 


{fhgi}(k„,r An ) = 0. 

(39) 


At the origin we impose lim^^o f/(k„, r An ) = 0 for the regular 
solution. At large distances the boundary condition of course 
depends on whether the energy e is positive or negative. For 
scattering neutron states (positive e), 

lim gi(k n ,r An ) e l(KrAr ~^\ (40) 

rAn~* 00 

while for final neutron bound states (negative s), 

lim 9l (k n , r n ) ^ (41) 

rAn-* 00 
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with k„ — \J-2m n e/h. 

If the imaginary part Wa„ of the neutron-target optical po¬ 
tential is small, we can use first order perturbation theory to 
express 


GT( r A , 


r A„; £ ) 


ft y <f > n( r A n ) , Pn( r An) 

2 m„ 2-J E n + iTJ2 - E 1 


(42) 


where 


r = 2 

1 n ^ 


f 


Anfyn d^Ani 


(43) 


and <p n ,E n are the eigenfunctions and eigenvalues of the 
Schrodinger equation corresponding to the real part of the op¬ 
tical potential. If we now consider an energy E close to an 
isolated resonance, i.e. \E-E„\ «: | E m —E„\, only the nth term 
of the sum will contribute to the Green’s function, and 


G° pt (r An ,r' A -k) 


ft K^An^niTAn) 
2m n E„ + IT nil - E' 


(44) 


The resulting neutron wave function is 


<A(r An) = 


ft <!>n( r An) 

2 m n E n + iT n / 2 - E 


f 


£(r' A n)S (r' An )dr' An . 


(45) 


According to the particular nature of the source term (see Eqs. 
( [23] ) and pT[), the integral in Eq. (45 i has the form of a one- 
neutron transfer DWBA amplitude 


't'(INT) 

1 n 


f f n S n dr' An = J f n {x ( p\v post , prior | hXi)dr' An 


(46) 


to the single-particle state <f>„ of the target-neutron residual 
nucleus. We can then write: 


j(1NT) 

ikniXAn) = - -- ~<Pn( r Ari), (47) 

E n + lT„/2 - E 

and the final neutron wavefunction il/„(r All ) can be inter¬ 
preted as the nth eigenstate of the neutron-target (real) single¬ 
particle potential times the direct transfer amplitude to this 
particular state, modulated by an energy denominator. The 
absorption cross section is proportional to the matrix element 


„ (1 NT) I ■ 


WWAnty) = 


(E„ - E) 2 + T 2 /4 

i r„ 


f* 


2 (E„ - E) 2 + T 2 /4 


W A n<Pn dr An 


|^ 1NT) | 2 


(48) 


As a consequence, if the transfer amplitude r, ( , 1NT) is approx¬ 
imately constant in an energy interval of the order of T„, 
the energy-dependent differential cross section around a reso¬ 
nance has a Lorentzian shape, and the integrated cross section 
under the peak is independent of W An for small enough WAn¬ 
il can be shown (see Appendix C) that there is a simple rela¬ 
tionship between the cross section for the capture of a neutron 
in a bound state of finite width and the cross section for direct 



V 

W 

W D 

a 

ClD 

r 

r D 

rc 

d 

99.0 

0.0 

16.7 

0.84 

0.64 

1.12 

1.31 

1.30 

P 

50.6 

0.0 

14.1 

0.678 

0.47 

1.25 

1.25 

1.25 


TABLE I: Optical model parameters. Energies are expressed in MeV 
and lengths in fm. In the proton channel, the parameters are those 
listed in m for proton- 93 Nb scattering at 16.2 MeV, while in the 
deuteron channel they correspond to those listed in the same refer¬ 
ence for deuteron- 93 Nb scattering at 17 MeV. 


transfer to the corresponding zero-width bound state. Assum¬ 
ing again that 7’/ ; INT) is essentially constant in an energy range 
of the order of T„ = 2(W A „), we have 


d 2 cr 

dQ.pdE p 


NEB 


(E, Q) 


1 I dcr n 

2 n (E n - E) 2 + I*/4 dn 


(O), (49) 


where is the direct transfer differential cross section to the 
«th eigenstate of the real potential. At the resonance energy 
peak (E = E n ), we have the simple relationship 


d 2 cr 

dYlpdEp 


(E = E n ,Sl) 


NEB 


2L<Ei (0) . 

r„7r dn J 


(50) 


III. RESULTS 
A. Numerical details 

We present the results obtained for the reaction 93 Nb(r/, p) 
at two different beam energies, E,i = 15 MeV, and = 25.5 
MeV. The optical model potentials U A d, U Ap used in the initial 
(deuteron) and final (proton) channels respectively were taken 
from the Perey and Perey compilation ( Ifl9l ), and are summa¬ 
rized in Table [T] For the final state interaction U An between 
the neutron and the 93 Nb target, we have slightly modified the 
energy dependent Koning-Delaroche global optical nucleon- 
nucleus potential ( IT8l ). by setting the real part parameter to 
V = 50.3 MeV in the whole energy range, and keeping the 
original energy dependence of the other parameters. This was 
needed to reproduce essential features of the 94 Nb nucleus. 
Furthermore, although we respect the energy dependence of 
the depth of the surface imaginary part Wf> of the Koning- 
Delaroche potential, we do not let it fall below 4 MeV, corre¬ 
sponding to the experimental energy resolution of the data in 
m. The maximum partial wave l p used in the calculations is 
15 and 20 for Ed = 15 MeV and Ed = 25.5 MeV respectively. 
The contribution of final neutron states with / > 8 is found to 
be very small. 

The deuteron ground state wavefunction is taken to be an 
L = 0 state with a radial wavefunction generated by a Woods- 
Saxon potential with radius Rd = 0.4 fm and diffusivity 
ad = 0.6 fm. When the real depth is adjusted to reproduce 
the binding energy of deuteron, the resulting wavefunction is 
compatible with the experimental value of the mean square ra¬ 
dius of the deuteron and the zero-range constant Z) () = -122.5 
MeV.fm 3/2 . 
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FIG. 2: Energy distributions of the detected proton for 93 Nb(d,p) 
at (a) 15 MeV and (b) 25.5 MeV: total cross sections (solid line), 
elastic breakup (EB, dot-dashed line) and inelastic breakup (NEB. 
dashed line). The arrows indicate the position of the neutron emis¬ 
sion threshold. Data is from ED. 


In Fig. [2] we show the proton energy distributions for 
93 Nb(d,p) at 15 MeV and 25.5 MeV along with the data fl6l . 
Also shown is the breakdown into elastic breakup (EB) and 
inelastic breakup (NEB). Our results indicate that the inelas¬ 
tic breakup is dominant at all energies, nevertheless elastic 
breakup is important particularly around the peak of the dis¬ 
tribution. The comparison with the data demonstrates that our 
model provides a good account of the process. A close com¬ 
parison with the theoretical results presented in M show sig¬ 
nificant differences that can be partially attributed to the ne¬ 
glect of the additional terms arising from non-orthogonality. 

We also compared our results with those from |[T3l . At the 
peak of our distribution for the higher beam energy ( E p = 14 
MeV), the elastic breakup contributes 25% of the total cross 
section. This differs significantly from the results presented 
in Q3), where the elastic breakup is less that 10% of the total 
cross section at the peak. This is an important issue because 
elastic breakup will not lead to the neutron being captured 
into a compound state and therefore will need to be subtracted 
from the total cross section in order to apply the surrogate 
method. In the work of oca, the elastic breakup is treated in 




FIG. 3: Energy distributions for 93 Nb(d,p) at (a) 15 MeV and (b) 25.5 
MeV: partial wave decomposition . 


a separate formalism, namely with the continuum discretized 
coupled channel method. In addition, given that the inelastic 
contribution is computed in the post formalism, the authors 
used continuum bins to address the convergence issues. A 
more detailed comparison between these two methods will be 
very useful. 

For the surrogate method, the spin distributions of the (d,p) 
cross sections are important in making the connection to neu¬ 
tron capture (see Section [T]). In Fig. [3] we provide the break¬ 
down in terms of the various angular momenta, for both beam 
energies considered. A strong peak is found around E p = 10 
MeV, for L - 3 corresponding to a narrow resonance in the 
neutron-target system. A much broader resonance is present 
in the L = 1 channel and therefore the peak in that compo¬ 
nent is less pronounced. For L — 0 and L - 2, we can see 
the signature of the neutron-target bound states, at high pro¬ 
ton energy. From the distributions in Fig. [3] it is obvious that 
the Weisskopf approximation, typically used in the analysis of 
surrogate reactions, is not valid. The cross section is strongly 
dependent on the angular momentum and closely connected 
to the internal structure of the composite final state. 

One attractive feature of the model here developed is the 
ability to provide, in a consistent framework, a description 
of neutron capture into the continuum and into bound states. 
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FIG. 4: Energy distributions for 93 Nb(d,p) at 15 MeV below the neu¬ 
tron threshold, for those partial waves that have bound states: (a) 
strength of the imaginary potential W = 0.5 MeV. and (b) strength of 
the imaginary potential W = 3.0 MeV. 



FIG. 5: Effect of the non-orthogonality and cross terms in the cross 
sections: (a) energy distribution for 93 Nb(d,p) at 25.5 MeV for an 
angle 6 = 10° and (b) angular distributions for both E n = -3 MeV 
and E„ = 5 MeV. 


The optical potential used to describe the n-target system is 
primarily based on fits to elastic scattering, appropriate for 
positive energies, but our results depend on the continuation 
of the potential to negative energies. In Figs. [2] and [3j for the 
highest proton energies, corresponding to the neutron below 
threshold, we keep Wq = 4 MeV, as stated above, to describe 
limited experimental resolution. In Fig. [4j by contrast, we 
explore the effects of smaller imaginary terms in Uah in the 
93 Nb(d,p) reaction, by changing this value to VV'/j = 0.5 MeV 
(a) and Wo = 3 MeV (b). For the smallest imaginary term, 
bound states appear as narrow peaks. Increasing the imagi¬ 
nary term increases their widths. An example of this can be 
seen for the L — 2 states: when the imaginary term is small, 
we can resolve the L — 2 spin-orbit partners, while as we in¬ 
crease the imaginary term, these two states become blurred to 
the point of becoming indistinguishable in Fig. [3] The arrows 
in Fig. [^correspond to the location of the bound states for the 
real potential. We find that that in addition to the spreading 
of the states, the introduction of the imaginary component can 
shift the peaks to higher energy. As discussed in section 
the limit of Wo = 0 MeV in our formulation corresponds 
to the standard DWBA transfer to bound states. We have in¬ 
deed used this fact to check our calculations by comparing our 


results in this limit to those produced by Fresco ll20ll . Com¬ 
plete agreement was found. 

The basis of the large controversy between the AV and the 
UT approaches stems from the non-orthogonality term, dis¬ 
cussed in detail in section IU bI It is therefore critical to under¬ 
stand the importance of this term in the calculations. In Fig. [5] 
we plot the inelastic breakup cross sections, with (solid line) 
and without (dashed line) the non-orthogonality term. The 
neutron energy distribution, for 6 p = 10°, is shown in panel 
(a) and the angular distribution, for two neutron energies, is 
shown in panel (b). The non-orthogonality term is most im¬ 
portant at the peak of the energy distribution but remains im¬ 
portant even when the neutron is captured into bound states. 
In terms of the angular distribution, the non-orthogonality 
term manifests itself mostly at forward angles. The results 
in Fig. [5] demonstrate the need for the inclusion of the non¬ 
orthogonality term in the prior-form formalism. The discrep¬ 
ancies found in the literature between the AV and the UT ap¬ 
proaches are, to a large extent, a consequence of the fact that 
UT neglects the non-orthogonality term. 
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IV. CONCLUSIONS 

We have derived, implemented and validated a practical 
theory for computing cross sections for inclusive processes 
of the type A(d,p)X, where only the proton is detected in the 
final state. This includes elastic breakup as well as all other 
inelastic processes. Our model is based on the assumption 
that the proton is a spectator and cannot excite the target. The 
same framework is able to describe these reactions across the 
full neutron energy range in a consistent manner. Our model 
is directly relevant to the application of the surrogate method 
for extracting neutron capture cross sections from (d.p) reac¬ 
tions. The present formalism can be generalized to describe 
the partial fusion of any cluster in a loosely bound projectile. 

We discuss in detail the post and the prior formalism. While 
the cross sections derived in the post-form are formally more 
elegant, they pose convergence challenges. This challenge is 
removed in the prior formalism. However, additional terms 
due to non-orthogonality arise. We show that these terms are 
important and cannot be neglected for all energies of the de¬ 
tected proton. We also predict the spin distribution of the re¬ 
sulting compound nucleus, an important ingredient for the ap¬ 
plication of the surrogate method to extract neutron capture 
reactions. We show that the cross sections are strongly de¬ 
pendent on the angular momentum, a factor that needs to be 
incorporated in the analysis of the surrogate data. Finally, we 
also make predictions regarding the relative contribution of 
the elastic breakup and all other inelastic processes. The elas¬ 
tic breakup process does not lead to neutron capture and needs 
to be subtracted from the total (d.p) cross section when relat¬ 
ing it back to (n,y). In this work, we have resolved a decades- 
old controversy, namely that between Austern el al. and Uda- 
gawa et al., promulgating respectively the post and the prior 
methods for describing the inclusive process A(d,p)X. 

While our method is able to describe the one example we 
have studied, it is important to perform a systematic study 
for various nuclei across the nuclear chart. Also, since our 
method makes consistent predictions for both the elastic and 
inelastic breakup, it would be very useful to have the mea¬ 
surement of elastic breakup for the same cases that have been 
studied with the inclusive experiments. 

In comparing our methods with other calculations, we 
found small differences with lfl6l . which can be partially at¬ 
tributed to the neglect of the non-orthogonality induced terms. 
More concerning, we found significant differences with the 
predictions of lfl3l . particularly in the elastic breakup cross 
section. Given that in El the elastic breakup was obtained 
through a different formalism, the comparison is not straight¬ 
forward. Nevertheless, a better understand on the source of 
the differences is desirable. Again, elastic breakup data for 
these cases would be useful. 

Obtaining a reliable cross section for processes of the type 
Afd.pjX is only the first step in the application of the sur¬ 
rogate method. The elastic breakup cross section needs to 
be subtracted from the total cross section, and the remaining 
cross section needs to be treated within a statistical approach 
to determine the fraction of neutrons that end up being cap¬ 
tured versus those that evaporate. The coupling to a statistical 


method is planned for the near future. 
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Appendix A: Discussion of EB and NEB in the context of flux 
conversation 


We rewrite the differential equation for the neutron wave- 
function eq. (|24| in the following form: 


(v 2 + U An -E)t P n ost =S post , (51) 


where 


S 


post ~ 


2 m„ c 

n 2 post ’ 


(52) 


and U An = ^prU An , E = ^ [e - E p - s A ). From ( J 51 } and 
its complex conjugate, and multiplying respectively by 
and if/ post we get (we will drop the post suffix) 


^,(V 2 + U A n-E)^ n = r n S 

if/ n (v 2 + 0 An - = 0-„5 + . (53) 


Substracting the above equations, and integrating over a large 
volume, we obtain 



Wn - ^y^n) dr n - 2i J' \lj/ n 

-J 


W /\/ 7 dv n — 

3(iJ/ n S f )dr n . 


(54) 


The first term can be cast into an outgoing elastic flux across 
the surface enclosing the volume, while the second term ac¬ 
counts for the non-elastic breakup. The above equation de¬ 
scribes how the flux generated by the right-hand side term is 
converted into an elastic and a non-elastic contribution. 


Appendix B: Partial wave decomposition of the differential cross 

section 


with the source terms having a similar decomposition 

S (r„; k p ) = ^ s lm i p (r n ; k p ) Y lm {6 Bn )Y lp , n (6), (59) 


ImL 


<Pimi p (r A „, kp) = J Gi(r An , r’ An ) s, mlp (r An \k p ) rf n dr' An 
1 / f' An 

= y I gi(k„, r An ) fi(k n , r' An )s hn , p (r' An \k p ) r' An dr' n 

+ft(k n ,r An ) I gi(k n ,r' An )si m i p (r' An ;k p )r' An dr' An . 

Jr An ) 

(60) 


Appendix C: Transfer cross sections to bound states 

In this appendix we show that the value of the integrated 
cross section under the peak in Eq.(49) is equal to the direct 
transfer cross section to a sharp bound state. Let’s consider 
the limit 


lim WWa*W> = 7t|7^ 1NT) |“ 6(E n - E), (61) 

Wa,-* 0 1 1 


and substitute (f 6 ~i~[) in Eq. (25 1 , 


d~cr 


dEl p dE p 


NEB 


m d m p k p 
4i r 3 fl 4 k d 


{tyn\ W An Wn) 


_ m d m p kp I (1MT)|2 _ 

~4n 2 tdk d Vn 1 1 " h 


(62) 


After partial wave decomposition, the multiple differential 
cross section can thus be written in the prior and post forms, 

d 2 cr 2n 


dClpdEp tw d 
where the post contributions are 


Pp(Ep) J] 


B' 


,pnor,post 
Im 


(55) 


B' 


,post _ 
lm 


= f I Z M v -J0p)\ 2 WAn(r A n) drp n , (56) 

J l 

l D 


and the prior contributions 

c or = f 

i P 

+ f \^<P^ p {rBn\kp)Y , \ l Wp)fw An (r A n) dr Bn 
J I 

1 D 


■23t 


fz 

^ i r 


Kl^kpWiJ (r Bn -,k p ) 


x Y l SJ0 p )Y l i n (O p )W An (r A „) dr Bn . (57) 

Here we have used the neutron partial wave functions, equiv¬ 
alently for prior or post, 

^„(r B „;k p)- '^<l>i m i p (r B ,i,kp)Y , m (6 Bn )Y l _ p m (6p)/r Bn , (58) 

l,m,lp 


where the subindices p,d refer to the proton and the deuteron 


respectively, and we have used the density of levels (12 1 . If 
we integrate over a vanishing interval 6E around E„ , we get 


dcr 

d£l„ 


NEB 


m d m p k p 
4 7T 2 tl 4 k,i 


|^ 1NT f 


(63) 


which is the DWBA transfer differential cross section. 

Moreover, in this limit the non-orthogonality term and the 
cross term vanish. Indeed, the functions ip™(r Bn -,kp), not be¬ 
ing affected by the propagator, do not exhibit any resonant 
behavior around the poles E n . Thus the second term in Eq. 
( |57j ) vanishes as W An —» 0. On the other hand, the third term 
behaves as 


lim - 

r„^o (E„ 


E) + iT n /2 


(64) 


when W An —> 0, and is equal to zero everywhere except in 
the zero-measure interval around E n = E, where it has the 
finite value -2 i. Its contribution thus vanishes after the en¬ 
ergy integration. In addition, the interaction U An used in the 
prior representation (see Section [HB] i coincides now with the 
real potential V A „, so the standard post-prior symmetry is fully 
recovered. 


10 













[1] J. E. Escher, J. T. Burke, F. S. Dietrich, N. D. Scielzo, I. J. 
Thompson, and W. Younes, Rev. Mod. Phys. 84, 353 (2012). 

[2] R. Elatarik, L. A. Bernstein, J. A. Cizewski, D. L. 
Bleuel, J. T. Burke, J. E. Escher, J. Gibelin, B. L. Gold¬ 
blum, A. M. Elatarik, S. R. Lesher, et al., Phys. Rev. C 
81, 011602 (2010), URL http://link.aps.org/doi/18. 
1183/PhysRevC.81.811682 

[3] A. Kerman and K. McVoy, Annals of Physics 122, 197 (1979), 
ISSN 0003-4916, URL http://www.sciencedirect.com/ 
science/article/pii/8883491679983814 

[4] T. Udagawa and T. Tamura, Phys. Rev. Lett. 45, 1311 
(1980), URL http://link.aps.org/doi/18.1183/ 
PhysRevLett.45.1311 

[5] N. Austern and C. M. Vincent, Phys. Rev. C 23, 1847 (1981), 
URL http://link.aps.org/doi/18.1183/PhysRevC. 
23.1847 

[6] X. H. Li, T. Udagawa, and T. Tamura, Phys. Rev. C 30, 
1895 (1984), URL http://link.aps.org/doi/18.1183/ 
PhysRevC.38.1895 

[7] T. Udagawa and T. Tamura, Phys. Rev. C 33, 494 (1986), URL 
http://link.aps.org/doi/18.1183/PhysRevC.33.494 

[8] M. Ichimura, N. Austern, and C. M. Vincent, Phys. Rev. 
C 32, 431 (1985), URL http://link.aps.org/doi/18. 
1183/PhysRevC.32.431 

[9] M. Ichimura, Phys. Rev. C 41, 834 (1990), URL http:// 


link.aps.org/doi/18.1183/PhysRevC.41.834 

[10] M. Hussein and K. McVoy, Nuclear Physics A 445, 124 (1985). 

[11] M. S. Hussein and R. Mastroleo, Nuclear physics A 491. 468 
(1989). 

[12] M. S. Hussein, T. Frederico, and R. Mastroleo, Nuclear physics 
A 511, 269 (1990). 

[13] J. Lei and A. Moro, Physical Review C (2015), submitted. 

[14] B. V. Carlson, R. Capote, and M. Sin, arXiv:1508.01466 [nucl- 
th] (2015). 

[15] T. Udagawa, Y. Lee, and T. Tamura, Physics Let¬ 
ters B ~ 196, 291 (1987), ISSN 0370-2693, URL 

http://www.sciencedirect.com/science/article/ 
pii/8378269387987337 

[16] R. C. Mastroleo, T. Udagawa, and M. G. Mustafa, Phys. Rev. 
C 42, 683 (1990), URL http://link.aps.org/doi/18. 
1183/PhysRevC.42.683 

[17] I. J. Thompson and F. M. Nunes, Nuclear Reactions for Astro¬ 
physics (Cambridge University Press, 2009). 

[18] A. J. Koning and J. Delaroche, Nuclear physics A 713, 231 
(2003). 

[19] C. Perey and F. G. Perey, Atomic data and nuclear data tables 
17, 1 (1976). 

[20] I. J. Thompson, Comput.Phys.Rept. 7, 167 (1988). 


11 



